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A field-theoretical description of the behavior of disordered, elastically isotropic, compressible sys- 
tems characterized by two order parameters at the bicritical and tetracritical points is presented. The 
description is performed in the two-loop approximation in three dimensions . The renormalization 
group equations are analyzed, and the fixed points corresponding to different types of multicritical 
behavior are determined. It is shown that the effect of elastic deformations causes a change in the 
regime of the tetracritical behavior of disordered systems because of the interaction of the order 
parameters through the deformation field. 
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It has been shown in [1] that the presence of frozen 
structural point defects gives rise to changes in the behav- 
ior of a system in both bicritical and tetracritical regions. 
According to the cited paper, the effect of (^-correlated 
impurities is essential only for Ising systems and leads to 
a decoupling of the order parameters at the multicritical 
points. Later [2], it was shown that elastic deformations 
cause a renormalization of the effective charges for the 
interaction of critical fluctuations, which lead to an in- 
crease in the interaction of order parameters and to a 
change in the type of multicritical behavior. Therefore, 
the study of the combined effect of elastic deformations 
and structural point defects in the multicritical region is 
of great interest. 

The purpose of this paper is to study the influence of 
striction effects on disordered systems whose phase dia- 
grams already contain multicritical points of bicritical or 
tetracritical type. In the first multicritical point 

is the point of intersection of two lines of second-order 
phase transitions and one line of first-order phase transi- 
tions, while in the second case, it corresponds to the in- 
tersection of four lines of second-order phase transitions. 
In the immediate vicinity of a multicritical point, the 
system exhibits a specific critical behavior characterized 
by the competition between the types of ordering. At a 
bicritical point, one parameter is displaced by another, 
whereas the tetracritical point allows a mixed phase with 
a coexistence of different types of ordering. Such systems 
[I] can be described by introducing two order parameters 
that belong to different irreducible representations. 

In structural phase transitions that occur in the ab- 
sence of piezoelectric effect, in the paraphase, elastic 
strains play the role of a secondary order parameter 
whose fluctuations are not critical in most cases [3, 4]. 
Since, in the critical region, the main contribution to 
the striction effects comes from the dependence of the 
exchange integral on distance, only elastically isotropic 
systems are considered below. 



The model Hamiltonian of a system has the form 
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Here, $(x) and ty(x) are the fluctuating order param- 
eters; itoi and i*02 are positive constants; t\ ~ \T — 
T c i\/T c i and ti ~ \T—T C 2\/T C 2, where T c \ and T C 2 are the 
phase transition temperatures for the first- and second- 

3 

order parameters, respectively; y(x) — u aa{x) where 

u a fj is the strain tensor; g\ and 172 are the quadratic stric- 
tion parameters; (3 is a constant characterizing the elas- 
tic properties of the crystal; and D is the space dimen- 
sion. In this Hamiltonian, the integration with respect 
to the components depending on the nonfluctuating vari- 
ables, which do not interact with the order parameters, 
is already performed along with the replica procedure of 
averaging over impurities. The properties of the initial 
system can be obtained in the limit m — > 0. The nonneg- 
ative constants 801, £02 and <5o3 describe the interaction 
of critical fluctuations via the impurity field. The inter- 
action of impurities with elastic deformations has linear 
character and, in the case of averaging over impurities, 
leads to the redetermination of the constants <5oi , S02 and 
£03 [5]. 

By going to the Fourier transforms of the variables in 
Eq. (1), we obtain the Hamiltonian of the system in the 
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form 
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Here, the components yo describing uniform strains are 
separated. According to [3], such a separation is neces- 
sary, because the nonuniform strains y q are responsible 
for the exchange of acoustic phonons and lead to long- 
range interactions, which are absent in the case of uni- 
form deformations. 

We define the effective Hamiltonian that depends only 
on the strongly fluctuating order parameters $ and \Pof 
the system as follows: 

exp{~Hl<S>,y]} = B J exp{-tf [<i>, *,</]} n rf % ( 3 ) 

If the experiment is performed at a constant volume, the 
quantity yo is a constant, and the integration in Eq. (3) 
is performed with respect to only the nonuniform strains, 
while the uniform strains do not contribute to the effec- 
tive Hamiltonian. In an experiment at a constant pres- 
sure, the term Pft is added to the Hamiltonian, with the 
volume being represented in terms of the strain tensor 



components in the form 

n = n [i + ^2u aa + ^2u aa u l3 p + o(u 3 )} (4) 

a = l 

and the integration in Eq. (3) is performed also with 
respect to the uniform strains. According to [6], the in- 
clusion of quadratic terms in Eq. (4) may be important 
at high pressures and for crystals with strong striction 
effects. As a result, we obtain 
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This Hamiltonian leads to a wide variety of multicritical 
points. As for incompressible systems, both tetracritical 
(v 3 + 12{ Zl z 2 -w lW2 )) 2 < (ui + 12(«?-«;?))(«2 + 12(z|- 

and bicritical 

(w 3 + 12(ziz 2 -wiw 2 )) 2 > («i + 12(z?-wf))(w 2 + 12(z|- 

behaviors are possible. In addition, the striction effects 
may give rise to multicritical points of higher orders. 

In the framework of the field-theoretical approach [7], 
the asymptotic critical behavior and the structure of the 
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phase diagram in the fluctuation region are determined 
by the CallanSymanzik renormalization-group equation 
for the vertex parts of the irreducible Greens functions. 
To calculate the (3 and 7-functions as functions involved 
in the CallanSymanzik equation for renormalized inter- 
action vertices tti, 112, u 3 , 6\, 82, 82, 51, 32, g^, 9^ or for 
complex vertices Z\, Z2, W\, W2, V\, v 2 , v 3 , 8\, 62, 6 3 , 
which arc more convenient for the determination of the 
multicritical behavior, the standard method based on the 
Feynman diagram technique and on the renormalization 
procedure was used [8]. As a result, the following expres- 
sions were obtained for the functions in the two-loop 
approximation: 
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It is well known that the perturbative series expansions 
are asymptotic, and the vertices of the interactions of 
the order parameter fluctuations in the fluctuation re- 
gion are sufficiently large for Eqs. (6) to be directly ap- 
plied. Therefore, to extract the necessary physical in- 
formation from the expressions derived above, the Padc- 
Borel method generalized to the multiparameter case was 
used. The corresponding direct and inverse Borel trans- 
formations have the form 

f(vi,v 2 , v 3 ,S 1 ,S 2 ,S 3 , z 1 ,z 2 ,w 1 ,w 2 ) (7) 

= c i ,.. il0 <^ 2 ^ 3 ^ 4 4 5 4 6 ^ 7 4 8 <^2 
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For an analytic continuation of the Borel transform of a 
function, a series in an auxiliary variable 9 is introduced: 
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and the [L/M] Pade approximation is applied to this se- 
ries at the point 9 = 1. This approach was proposed and 
tested in [9] for describing the critical behavior of sys- 
tems characterized by several vertices corresponding to 
the interaction of the order parameter fluctuations. The 
property [9] that the system retains its symmetry under 
the Pade approximants in the variable 9 is essential in 
the description of multivertex models. 

In the two-loop approximation, the functions were cal- 
culated using the [2/1] approximant. The character of 
the critical behavior is determined by the existence of a 
stable fixed point satisfying the set of equations 

/3i(v*,V2,vl,S*,S2,S3,zl,Z2,wl,W2) = 

(i = l,...,10). (8) 

The requirement that the fixed point be stable is reduced 
to the condition that the eigenvalues 6, of the matrix 

d(3i{v$, v%, v$, 51, 6%, 5%, zl,z^wl,wl) 

= 5 

dVj 

(Vi,Vj = V*,V2,V2,5{,52,53,Zi,Z2,W*,W2) (9) 

lie on the right-hand complex half-plane. 

The resulting set of summed [3 functions contains a 
wide variety of fixed points lying in the physical region 
of the vertex values with Vi > 0. 

A complete analysis of the fixed points, each of which 
corresponds to the critical behavior of a single order 
parameter, was presented in our recent publication [5]. 
Now, we consider the combined critical behavior of two 
order parameters. 



The analysis of the values and stability of the fixed 
points offers a number of conclusions. The tetracritical 
fixed point of a disordered incompressible system (vi = 
v 2 = 1.58892, i; 3 = 0, tfi = 5 2 = 0.03448, 5 3 = 0, z x = 0, 
z 2 = 0, w\ = 0, w 2 = 0) is unstable under the effect 
of uniform deformations (b\ = b 2 = 0.461, 63 = 0.036, 
b 4 = b 5 = 0.461, b 6 = 0.036, b 7 = b$ = b 9 = b w = 
—0.236). The striction effects lead to the stabilization of 
the tetracritical fixed point of a compressible disordered 
system ( Vl = v 2 = 1.58892, v 3 = 0,5t = 5 2 = 0.03448, 
5 3 = 0, Z! = 0.04599, z 2 = 0.568836, w x = 0.017759, 
w 2 = 00.551849, h = b 2 = 0.461, b 3 = 0.036, b 4 = b 5 = 
0.461, b 6 = 0.036, 67 = 1-189, 6 8 = 0.003, 69 = 5.391, 
67 = 0.999). 

The question about the stability of other multicriti- 
cal points cannot be resolved in terms of the described 
model, because the calculations lead to a degenerate set 
of equations. The degeneracy can be removed by con- 
sidering the Hamiltonian with allowance for the terms of 
higher orders in both the strain tensor components and 
the fluctuating order parameters. 

Thus, the striction-induced interaction of the fluctuat- 
ing order parameters with elastic deformations, as well 
as the introduction of frozen point impurities into the 
system, leads to a change from the bicritical behavior 
to the tetracritical one. Elastic deformations lead to a 
change in the regime of tetracritical behavior of a dis- 
ordered system because of the interaction between the 
order parameters through the acoustic phonon exchange. 
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